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Continuous center-of-mass position measurements performed on an interacting harmonically 
trapped Bose-gas are considered. Using both semi-analytical mean-field approach and completely 
quantum numerical technique based on positive P-representation, it is demonstrated that the atomic 
derealization due to the measurement back action is smaller for a strongly interacting gas. The nu- 
merically calculated second-order correlation functions demonstrate appearance of atomic bunching 
as a result of the center-of-mass measurement. Though being rather small the bunching is present 
also for strongly interacting gas which is in contrast with the case of unperturbed gas. The per- 
formed analysis allows to speculate that for relatively strong interactions the size of atomic cloud 
determined with a single snapshot measurement can become smaller than the ground-state cloud 
size. 
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I. INTRODUCTION 

Now that generation of the trapped degenerate gases 
has become almost routine procedure in many labora- 
tories the experimental and theoretical activities are di- 
rected towards applications of these systems. New fre- 
quency standards [l| and various types of weak-force de- 
tectors @, H[ belong to such applications. Another im- 
portant application of such quantum systems is testing 
quantum-mechanical predictions in mesoscopic regimes. 

All these experiments involve measurements of some 
property of the trapped gas. As is extensively discussed 
in quantum theory, the measurements of small objects 
do inevitably disturb them even if they are prepared as 
so called quantum non-demolition [i[ ones. This issue 
becomes especially important if the measurement is con- 
tinuously performed during an experiment. Recently a 
few experiments have been reported where the condi- 
tions required for continuous BEC measurement have 
been realized 5-8]. One of these works, Ref. Q, es- 
pecially addresses the measurement back-action on the 
trapped non-interacting gas. Detailed understanding of 
measurement back-action is very important for success- 
ful realization of quantum-limited feedback control of the 
trapped gas [H, [l(| . 

An important degree of freedom is the collective or 
center-of-mass position of the trapped gas. This degree of 
freedom can be accessed via bringing the system into in- 
teraction with a few-mode external field. Hence, this type 
of measurements is conceptually simpler than a measure- 
ment resolving internal structure of the trapped cloud. 
Moreover, the collective degree of freedom is proven to be 
experimentally accessible in recent works cited above 0- 
8] . These facts motivate the theoretical analysis of quan- 
tum motion of a degenerate gas with the continuously 
measured center-of-mass position that we present in this 
article. 



Unlike some other works @, , where similar prob- 
lems are considered, we focus on a gas of interacting par- 
ticles. It is known that the center-of-mass (CM) motion 
of a harmonically trapp ed g as is not coupled to the rela- 
tive degrees of freedom [lj|. This means that the inter- 
particle interactions do not affect the quantum dynamics 
of the center of mass. The CM measurement does, how- 
ever, influence the properties that depend on the relative 
motion such as particle density distribution [Io| . These 
properties are also affected by the inter-particle interac- 
tion. Thus an interesting question arises: does the si- 
multaneous action of the CM position measurement and 
the interaction result in some nontrivial dynamics of the 
gas? In other words, how does the interaction strength 
influence the behavior of the measured gas? It should be 
mentioned that interesting effects in a different system 
that also contains an open interacting BEC have recently 
been discussed in Ref. [13 ] . 

Considering a harmonically trapped repulsively inter- 
acting ID Bose gas we show that interaction-induced 
nonlinearity provides a mechanism that partially com- 
pensates measurement back action and stabilizes the 
cloud spreading. This can be accompanied by the nar- 
rowing of the instant density profile of the trapped gas 
compared with the ground-state density profile. Note 
that the master equation used in this article to describe 
the effect of the measurement also describes situations 
where the CM of the gas is weakly coupled to a hot reser- 
voir [la ]. To some extent such coupling might describe 
quantum fluctuations of the trapping potential that are 
always present in an experiment. Thus the discussed ef- 
fects might be relevant even for the cases where there is 
no explicit observation or control. 

We base our discussions on the numerical solution 
of the many-body quantum problem using positive P- 
representation of the density operator [l4| • Although this 
approach is known to be limited to relatively small evo- 
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lution times |l5[ it seems to be least resource consuming 
allowing to obtain reasonable results on a single desk- 
top computer. As a complement to numerical results we 
derive simplified mean-field approach that qualitatively 
agrees with our numerical results. 

The article is organized as follows. In Sec.|IT]the model 
that is considered is described and examples of possible 
physical implementations of the CM position measure- 
ments are given. The mean-field approach is derived and 
discussed in Sec. Mil Section |IV] is devoted to the nu- 
merical analysis of the system and to the comparison of 
the results with the mean- field predictions. In Sec. [V] we 
present the summary of the results of the article. 



II. MODEL 

Let us consider a gas of spin-0 bosons with mass m 
interacting via a repulsive ^-function potential. The gas 
is trapped in a cylindrically symmetric potential along 
x axis with frequencies ujq and lj±_ in axial and ra- 
dial directions, respectively; ao = (mcjo) -1 ^ 2 and = 
(muj^)^ 1 / 2 [25| are the corresponding lengths of the har- 
monic oscillator ground state. In the case of tight con- 
finement in radial direction (ujq <C and at low enough 
temperature all atoms are in the ground state of the ra- 
dial potential and the gas can be considered as being 
effectively a trapped ID gas. In quasi-lD case given that 
a_L ^ &0 the interaction of bosons can be described by 
the coupling constant giD — 2a s uj± [IH, where a s is the 
s-wave scattering length. 

We assume that the center-of-mass or, more precisely, 
collective position of the gas in a trap is continuously ob- 
served. Strictly speaking one should distinguish between 
the center-of-mass and collective position [2| . The former 
is the true mean value that accounts for particle number 
fluctuations while the latter is merely the weighted sum 
of particles coordinates. Below we will deal with states 
having small particle number fluctuations and neglect the 
mentioned difference. 

As a physical model of such a measurement one can 
suggest, for example, the approach discussed in [l7[ and 
experimentally realized in In these experiment the 

87 Rb atoms trapped in an atom chip have been coupled 
to high-finesse optical cavity. The authors show that the 
atoms can be localized entirely within a single antinodc 
of the standing- wave cavity field. Under these conditions 
the collective atomic position determines the atom-field 
coupling strength and can be determined by measuring 
the cavity transmission. 

Another straightforward example is Faraday rotation 
of the polarization of non-resonant light passing through 
the gas placed in a non- homogenous ma gnet ic field. 
Within the framework of the classical theory [19[ one can 
easily find the following expression for the total rotation 
angle 
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where e is the charge of an electron, m e is the mass of 
an electron, c is the velocity of light in a vacuum, uj is 
the probe light frequency, wr is the atomic resonance 
frequency, n(x) is the atom density, B(x) is the mag- 
netic field. Here we integrate over the atomic sample 
thickness and neglect a deviation of the refractive in- 
dex from unity. It is seen from Eq. ([TJ that in the case 
of the linear magnetic field the Faraday rotation is pro- 
portional to the center-of-mass coordinate of the atoms 
X = (1/-/V) J dxxn(x), where N is the estimate of the 
number of atoms that, for simplicity, might be thought 
of as the average number of atoms. One certainly can 
think of other schemes to measure the CM coordinate. 
Therefore, below we do not concentrate on a particular 
realization but consider the effects typical for all possible 
experimental approaches to CM position measurements. 

Our main concern is the measurement back action, 
thus we ignore particular measurement outcomes and an- 
alyze unconditioned dynamics. In this case the system is 
described by the master equation 



-i[H,p]- K [X,[X,p]]. 



(2) 



The first term on the right-hand side of this equation cor- 
responds to the hamiltonian evolution. In second quan- 
tization the Hamiltonian has the form 



H 



1 



J 2 



2m 

g m &(x) 2 ${xf]dx, 



l> t (x)<i>(a 



(3) 



where <t(x) and $t (x) are bosonic field operators obeying 
the commutation relation <&'(x')] =S(x — x'). The 

second term in Eq. ^ describes the measurement of the 
CM coordinate of the atoms 



X = ^ J x& {x)®{x)dx. 



(4) 



n(x)B(x)dx, 



(1) 



The parameter k characterizes the measurement resolu- 
tion of the apparatus and determines the back-action 
strength. Note that the same master equation can be 
obtained if all the atoms are weakly coupled to the same 
heat bath of high temperature. Namely, taking these 
limits into account Eq. ((2J follows from the well known 
Caldera-Legget master equation [la ]. 



III. MEAN-FIELD APPROXIMATION 

The evolution of the continuously measured interact- 
ing Bose gas can be found from the Hamiltonian ([3]) and 
the master equation ([2]). However, even for the moder- 
ate numbers of atoms N the direct numerical integration 
of the master equation is impracticable due to the large 
dimensionality of the N-atom Hilbcrt space. This prin- 
ciple problem of many-particle physics can be attacked 
with various numerical methods. In this article we use 
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approach based on positive P-representation of the den- 
sity operator However, before presenting numerical 
results we give a simple mean-field consideration that 
qualitatively describes the system behavior. It is shown 
that even this approach reveals some interesting features 
in the dynamics of the continuously measured interacting 
Bose gas. 

We start by defining a single-atom density matrix as 



(5) 



The evolution of the single-atom density matrix is de- 
scribed by the following equation 



p( Xl ,x 2 ) = \i^(dl 2 -dl x ) 



(x\ x\) 



+ 2ig 1D [n(x 1 ) - n(x 2 )} 

~ j^2~( X2 ~~ x 1 ) 2 ^p{x 1 ,x 2 ), 



(6) 



which was obtained from the master equation ([2]) ne- 
glecting density-density correlations. More precisely, the 
following approximation has been used 

(^( Xl ) 2 ^( Xl Mx 2 )) « «n(asi)) - ^-)p{xx,x 2 ). (7) 

Note that in the absence of the measurement, that is 
without the last term in the right-hand side of Eq. ([5]), 
this equation is a simple generalization of the well- 
known Gross-Pitaevskii equation (GPE) for the conden- 
sate wave-function. This is easily seen by substituting 
the coherent state p(x\,x 2 ) = (p*(xi)ip(x 2 ) into Eq. 
The measurement, as can be seen from Eq. ©, results in 
decay of non-diagonal elements (coherence) of the single- 
atom density matrix [2(| ■ This measurement- induced de- 
coherence prevents us from describing the gas in terms 
of the condensate pure-state wave-function. 

Using Eq. ^ the evolution of single-atom fluctuations 
can be derived as 



t (Ax 2 ) = -({x,p}), 
m 

d t (Ap 2 ) = -mu 2 ({x,p}) - 4g 1D (n'(x)p) + ^ 
9t{{x,p}) 



(8) 



2k 



N 2 ' 



— (Ap 2 ) ~ 2muj 2 (Ax 2 ) - 4g 1D (n' (x)x) , 
m 

where ({x,p}) denotes the anticommutator of x and p, 
n'(x) denotes the derivative of the density distribution 
with respect to the coordinate. Here we restrict the con- 
sideration to the states with (x) = (p)=0. The system (O 
is not closed since it contains the terms proportional to 
(n'(x)p) and (n'(x)x) that in general cannot be expressed 
via the single-atom second moments only. 

To render the system ([8]) closed we perform the follow- 
ing rough approximation, which, as will be seen below, 
is enough to grasp its qualitative behavior. First we note 
that without interactions the continuous measurement of 
the CM coordinate does not change the shape of the dis- 
tribution function describing the state of the atoms. This 



follows directly from the solution of the Fokker-Planck 
equation (FPE) for the Wigner function W(x,p) in the 
case of the non-interacting gas 



d t W(x,p) 



o x 



—a 2 

N 2 p 



W(x,p). (9) 



Solving this equation analytically with Gaussian initial 
condition one easily finds that the measurement of BEC 
only changes the width of the distribution preserving its 
Gaussian shape. We assume that in the considered situ- 
ation the distribution also remains approximately Gaus- 
sian during the system evolution. This is certainly not 
true for the strong interaction case [ll|. However, for 
weak interactions this might be a rather good approxi- 
mation. In this case one obtains the following expressions 
for the sought averages 



(n'(x)p) 
{n'(x)x} 



N 



8y/^(Ax 2 ) 3 / 2 
N 

Ay/^iAx 2 ) 3 / 2 
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(10) 



Substituting this result into the system ([8]) one obtains 
the following closed system of equations 



d t (Ax 2 ) 
d t (Ap 2 ) 



~{{x,v}), 
n 

-mil 2 cS ({x,p}) + 
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d t ({x,p}) = -(Ap 2 )-4mr! c 2 ff (Ax 2 ), (11) 



where the effective frequency fi e ff defined via 
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2v^m(Ax 2 ) 3 / 2 



(12) 



has been introduced. The effective frequency is deter- 
mined by the size of the atomic localization domain, 
which is represented by (Ax 2 ). The set of nonlinear ordi- 
nary equations (fTT|) can easily be solved using one of well 
established numerical techniques [2l[ . We use predictor- 
corrector Adams scheme that is known to be well suited 
for stiff problems. 

To characterize the effect of the CM position mea- 
surement on the gas we introduce the so called relative 
spreading of the atoms rj, defined as 



' no-mcas 



(13) 



Here, the subscripts "meas" and "no-meas" are used to 
distinguish the cases with and without the measurement, 
respectively. Taking for the initial values of the fluctua- 
tions the results obtained from the solution of the time- 
independent GPE we obtain the results shown in Fig. [T] 
In this plot the dependence of the relative spreading on 
the interaction parameter gir>N is shown for two values 
of measurement strength k/N 2 . Hereinafter, length is 
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FIG. 1: Relative spreading of the atomic cloud (rj) as a func- 
tion of the interaction strength for k/N 2 = 1 (blue squares) 
and k/N 2 =5 (red circles). 



measured in units of ao, time in units of Tq/2tt, where 
To is the period of oscillations in the harmonic potential. 
Thus, gir>N is the dimensionless quantity. The results of 
Fig. [T] correspond to the time instant t = To/ 4, which is 
the time when initial momentum uncertainty transforms 
to the position uncertainty. It is seen that being always 
positive the relative spreading r\ goes down for larger val- 
ues of the interaction constant. This means that the CM 
measurement always increases the width of the atomic 
localization domain, but for stronger atomic interactions 
the measurement-induced spreading is slower. Thus, an 
interacting trapped gas appears to be more stable against 
the measurement back-action noise than an ideal one, at 
least in the beginning of the noise-governed evolution. 

A simple explanation of this effect might be obtained 
on the basis of Eq. (fl~2|) . The effective frequency repre- 
sents a degree of atomic localization in an effective po- 
tential that is a combination of the trap and mean-field 
potentials. According to Eq. (TT21 the measurement in- 
duced spreading of the cloud due to non-linear response 
also increases the effective frequency, which corresponds 
to better localization. This mechanism partially compen- 
sates for atomic derealization due to the measurement 
back action. As follows from Eq. (fT2")) the effect should be 
more pronounced for larger atom-atom coupling. How- 
ever, the validity of this approach is certainly limited by 
5id N < 2 0Fmwg (Ax 2 ) 3 / 2 . 

There is another interesting feature of the dynamics 
of an interacting gas subjected to the CM position mea- 
surement. Note that the CM of harmonically trapped 
gas is not coupled to the relative motion [llj. Thus, 
the measurement-induced derealization of the CM is in- 
creased regardless of atomic interaction strength. The 
single atom distribution on the contrary depends on the 
interaction and for strong interaction the measurement- 
induced change of this quantity can be rather small. 
This implies appearance of correlations corresponding to 




FIG. 2: Two-atom illustration of the measurement-induced 
narrowing of the instant atomic density profile. The oval sym- 
bolizes the region of a single-atom localization that remains 
almost unchanged. The pulsed-shape distribution represents 
the CM position uncertainty that gradually increases due to 
the measurement back action. These two conditions can be 
simultaneously fulfilled if the atoms tend to bunch together. 



bunching of atoms and narrowing of the instant density 
profile compared with the initial unperturbed one. 

The easiest way to demonstrate this is to consider only 
two atoms. Assume that initially each of them spread 
over the same region in the trap (Fig. [2]). For the sake 
of simplicity let the atoms be spatially anti-correlated 
(perform anti-phase oscillations). In this case, the size of 
the atomic cloud is equal to the size of the single-atom 
localization region, while the CM is strongly localized 
near the center of the trap. The unconditioned mea- 
surement of the CM leads to its derealization, whereas 
the size of the atomic cloud, as predicted by the the- 
ory described above, remains almost unchanged. As il- 
lustrated in Fig. [5] to ensure such dynamics the atoms 
have to get closer to each other. This might result in 
an instantaneous density profile that is narrower than 
the initial ground-state one. In an experiment this effect 
manifests itself in a size of a single (not ensemble aver- 
aged) resonance-image or contrast-image snapshot that 
should become narrower after some time of the system 
evolution. 



IV. NUMERICAL SIMULATIONS 

In this section we discuss the results of ab initio nu- 
merical simulations of the system. We use the numerical 
scheme based on the positive P-representation [TJ. In 
spite of known drawbacks [l5| of this approach it is ef- 
ficient for dynamical calculations restricted to relatively 
short evolution times or small nonlinearities. In addition, 
this method is relatively easy to implement. 



A. Phase-space representation 

To treat the problem numerically the continuous 
atomic distribution is approximated by a lattice model. 
The space region occupied by the atoms is divided into 
M equal cells of length Ax: Xi = iAx, where i = 1, . . . , M. 
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For each cell i we define annihilation and creation oper- 
ators 



Xi+Ax/2 

i-Ax/2 
Xi+Ax/2 

Xi-Ax/2 



dx&(x)/VA~x, 



dx$ t (x) I \TKx , 



(14) 



that obey commutation relations [04, a,j] = and [&£, at] = 
Sij. Then the Hamiltonian ((3J) can be approximated by 
the following Bose-Hubbard Hamiltonian 

h = Y j ^A- ] + ^YM) 2 ^ ( 15 ) 



Here the matrix T accounts for the kinetic energy and 
potential energy in the trap, it is defined as 



1 



Ax 



x 2 \ 1 1 

' ~^l §ij -2A^ 6i+lj -2A^ 6i - lj - (16) 



The discrete version of the non-Hamiltonian term of the 
master equation @ is obtained in the similar way and 
need not be explicitly written. 

The total (many-atom) density operator is expanded 
using a positive P-rcprcscntation 

P (+) 



p(+>(a)A(a)d 2A W M /3 



using the following operator basis 

|a></3* 



A(a) = 



</3» ' 



(17) 



(18) 



Here, a = (a, /3), where a = {ax, . . . , cxm } an d (3 = 
{/?!,..., /?m} are complex vectors with components on = 
a'i+ia" and /?$ = /3^ + i/3". The P( + ) representation is 
guaranteed to always produce positive-definite diffusion, 
which is a necessary requirement for a stochastic differ- 
ential equation. 

Substituting the expansion <\17\ in the master equa- 
tion an d using the standard operator identities 



aA(a) = aA(a), 

atA(a) = (0 a +/3)A(a), 

A(a)a = (dp+a)h(a), 

A(a)a f - /3A(a) 



(19) 



one can obtain the following FPE for the positive PW 
function 



(9/ 



P<+>(a) = 



- + ^^(BfartBLW (20) 



B 



U^P (+) (a). 



Here, <9 M denotes partial derivatives d/dot^ if ju < M 
and dfdfi^-M otherwise, /i and ^ take values = 



1,...,2M. The elements of the drift vector A 
{Ax, A M+l ,. . .} are given by 



A, 



■it -9lD 2o K 2 im\ 

-tTyay - l ^T"i Pi " J^2 X i a ^ ( 21 ) 

3lD „2 



Ai+M = iTj^j + a» - Tn x i& 



N 2 ' 



The diffusion matrix can be divided into two parts cor- 
responding to different noise sources acting on the atoms. 
One of these sources may be attributed to the atom-atom 
interactions. It is described by the diagonal matrix 



B 



int 



9lD 



diag{(l -*)«!,..., (l + i)/3i,...}. (22) 



The other noise source is the measurement back action. 
This noise is represented by the matrix B moas with only 
one (first) non-zero column. This matrix is written as 



B, 



/2k 



(x\OL\ 



V 
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(23) 



It is easy to check that the total diffusion matrix is given 
by D = B int B mt + B moas Bm 0as - Assuming that the noise 
sources discussed above are represented by independent 
Wiener processes one can show that the FPE (|2"U1) is 
equivalent to the set of Ito stochastic differential equa- 
tions 



da = A{a, t)dt + B int {a,t)dW int {t) 



B, 



(a.,t)dW meas (t). 



(24) 



For the numerical simulations it is instructive to obtain 
a set of equations for real functions instead of Eq. (|24p . 
To do so we following, for example, Ref. [22j decom- 
pose the drift vector and the noise matrices into real and 
imaginary parts as A = A' + iA" and B = B' + zB". 
Since A(a) is an analytic function the derivatives d a 
and dp can be evaluated in either real or imaginary di- 
rections so that the resulting drift and diffusion terms 
can always be made real. Taking this into account 
one can define the new 4M-dimensional real drift vec- 
tor A = {A[, A' M+1 , . . . , A'{, . . . , A'm +1 , . . .} with the 
elements 



A' 



T ^ + if K ' a « + n ^ ,} - iv^' (25) 



A'! = -T 



9m 



(II If f/\ Z ft 

— (n^ - n t a 4 ) - — ^ , 



A' 



i+AI 



A" 



Ax 



N 2 

K 



= T4 + &iA'-n5W)-^ 



Ax 



N 2 



Here, r^ = a^-</3f and < = «•/?■' +a"/3- are real and 
imaginary parts of the complex atom number m = n'^Hn" . 



The new stochastic matrices B int and B mcas are 



B[ n 
B" 



(26) 
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with 



B. Results of numerical simulation 



B',, 



/SID,. r , 

A^ diag{ai 



B int 




ll 



+ #',... },(27) 



and 



with 



B, 



n; r: 

B" 



9/ 



B 



/2k 



iV 



and 



V 



/ Xia^ 



b: 



V 



(28) 
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(29) 



(30) 



The matrix in Eqs. flUJ and (|28j denotes the 2M x 2M 
matrix with zero elements. 

The SDE (j2~4")l is then cast into the following form con- 
taining new real 4M-dimensional Wiener noise vectors 



W int and W meas : 



da = A(a,t)dt + B iat (a,t)dW iat (t) 



B, 



3 (M)rfW meafl (t). 



(31) 



The 4M dimensional real vector a is formed of real and 
imaginary parts of a and (3. The elements of noise vec- 

with the elements denoted by 
obey the following properties 



tors dW int and dW meas 
dW, (int) and «, (meas> 



The numerical solution of SDE (|3"Tj) has been found 
using semi- implicit method discussed in Ref. 23]. This 
approach is generally quite stable when applied to non- 
linear and/or stiff problems. However the problem under 
consideration seems to have intrinsic instability which 
manifests itself during the simulations regardless of the 
stability of used numerical technique. Such an instabil- 
ity, as discussed in Ref. [Hj], is a result of incorrectly 
ignored boundary term during the derivation of FPE for 
the positive P-function. 

The conditions when the problem becomes unstable 
can be grasped looking at the terms of SDE pro- 
portional to the interaction constant. As soon as the 
system evolves to a quantum state with a* ^ (3 these 
terms may acquire a real part responsible for exponen- 
tial growth. Thus considerable (in some sense) deviation 
of the (3 from a* indicates the limits of applicability of 
the approach. 

In a series of numerical experiments performed for 
different values of the interaction constant the evolu- 
tion time has been determined during that the dynamics 
demonstrates no sign of "exploding" trajectories. This 
"secure" time interval is about quarter of the trap os- 
cillation period. Thus the numerical simulation results 
presented below correspond to To/4. The spatial dis- 
cretization used for the calculation is Ax = 0.33, the time 
step is Ai = 10~ 4 . Further decrease of these values does 
not change the appearance of the plots presented below. 
Thus the discretization error can be estimated as being 
of the order of the line thickness. The number of the 
stochastic trajectories equals to 20000. The estimated 
sampling error in this case is about 0.03% which is also 
below the thickness of the plot lines. 

For the initial state of the system it is convenient to 
take BEC broken-symmetry coherent state. This state 
for the lattice model reads 



I*) 



(33) 



(dW 



(mcas) \ 



(<W. 



(int) dw (int) ) 



0, 

Sijdt, 



(dWi dW) ') = Sijdt, 



(dwf nt W mcas) 



0, V ij. 



(32) 



where ip Xi is the value of the solution of the time- 
independent GPE in the space point Xi. The state gives 
the following initial values for the phase-space variables 



a'M = Re(^), #(0) = a((0), 
oj'(0) = Im(^ Xi ), A"(0) = oT(0). 



(34) 



Note that the measurement noise matrix B mcas consists 
of a single non-zero column. This means that all the 
modes (lattice cells) of the trapped gas are affected by 
the same measurement-induced noise. This is expectable 
since the considered noise acts on the collective (CM) 
degree of freedom of the system. The noise matrix B int 
that originates from the atomic interactions is diagonal, 
which means that each space point of the gas is driven 
by its individual Wiener noise. Noise sources acting on 
different coordinates of the gas are statistically indepen- 
dent. 



The results of numerical calculations of the atom den- 
sity profile for different values of the interaction strength 
giiyN are shown in Fig. [3] The curves shown by empty 
and filled symbols correspond to the situations with 
and without the measurement, respectively. For small 
interaction strength the density profiles corresponding 
to these cases differ substantially (circles). The differ- 
ence becomes less essential with the increase of guyN 
(squares). For yet stronger interaction the effect of 
the measurement becomes practically negligible, compare 
curves plotted with filled and empty triangles. 
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FIG. 3: Atom density profiles for different coupling constants. 
The curves shown by empty and filled symbols correspond 
to the situations with and without the measurement, respec- 
tively. 



Figure @] compares the relative spreading of the atoms 
V Q13|) as a function of the interaction strength for nu- 
merical simulation (filled symbols) and the mean-field 
approach of Sec. IIIII (empty symbols). Curves shown by 
squares (blue) and circles (red) correspond to different 
values of k/N 2 . It is seen that both methods predict 
similar qualitative behavior of the relative atomic spread- 
ing. That is the numerical simulations confirm predicted 
earlier decrease of spreading rj with increasing of atom- 
atom interaction strength. However, some quantitative 
discrepancy is observed, which is found to be more pro- 
nounced for larger gijyN. This is not surprising since 
in deriving the mean-field approach a couple of not- well- 
justified assumptions have been made. One of them is 
the assumption that the density-density type correlations 
are small. This is exactly the case for the chosen initial 
state ([55]) but may become wrong after some time of the 
system evolution. The other assumption that is certainly 
violated for a strongly interacting gas is the gaussian den- 
sity profile approximation. These poor approximations of 
the mean-field approach result in quantitative difference 
between the results of the two methods. Nevertheless, the 
essential feature of interacting gas dynamics subjected to 
the CM position measurement can be grasped within the 
mean-field approximation as derived above. 

The qualitative arguments presented at the end of 
Sec. |TTT] indicate the possibility of generating bunched 
states. These states are also characterized by squeezed 
compared with the unperturbed BEC ground state den- 
sity profile. Clearly the discussed bunching effect should 
manifest itself in the second-order correlation function, 
which can be defined as 

($t(0)$t(,x)d(.x)$(0)) 

g2\x) = — . loo) 

(*t(0)*(0))($t(a;)*(x)) 

This quantity characterizes density-density correlations 
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FIG. 4: Relative spreading of the cloud {if) as a function 
of the interaction strength. Filled and empty symbols show 
the results of the numerical simulation and the mean-field 
approximation, respectively. Squares (blue) and circles (red) 
correspond to different values of the measurement resolution. 
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FIG. 5: Second-order correlation functions for k/N 2 — 5 
and various values of the atom-atom interaction strength 
gmN. 



between the trap center and the point with the coordi- 
nate x. In case of independent densities in these points 
gi{x) — I, increased (decreased) likelihood to detect 
atoms separated by x means g% (x) > 1 (g%{x) < I). 

Figure [5] shows these density-density correlations for 
k/N 2 — 5 and different values of atom-atom interactions. 
It is seen that for small coupling strength giuN the mea- 
surement of the collective coordinate leads to the well 
observed bunching of the atoms (the curve shown with 
circles). This indicates that inside the space region oc- 
cupied by the atoms they are distributed non-uniformly, 
that is the atoms are grouped in a bunch with the size 
smaller than that of the occupation region. This is a 
natural result taking into account that the measurement 
back-action noise acts on the collective observable result- 
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ing in the overall delocalization of the atoms while keep- 
ing the initial instant density distribution. 

For larger coupling strength, the curve shown with 
squares in Fig. [5j the value of the second-order corre- 
lation function g2{x) is smaller for the small separations, 
x — > 0, than in the previous case. This indicates the 
smaller degree of bunching that is the result of increased 
atom-atom repulsion. For yet stronger repulsive inter- 
action (the curve shown with triangles in Fig. [5]) the 
value of the second-order correlation function 52(0) is 
only slightly larger than one (cannot be seen on the plot). 
Thus, the atom density in different locations becomes in- 
dependent. This is in contrast with the closed strongly 
interacting gas which exhibits anti-bunching [24|. 

It is worth noting that for weak (blue curve) and mod- 
erate (red curve) interactions the second-order correla- 
tion function does not seem to approach unity even for 
large separations. We believe that this is the result of col- 
lective character of considered back action. That means 
that essentially all the atoms are simultaneously sub- 
jected to the same back-action noise, which results in 
the observed in Fig. [5] long-range correlations. 

However we consider a finite system where the limit of 
very large separations is somewhat ambiguous. At least 
in numerical calculations of the second-order correlation 
function we have to restrict the separations to be smaller 
than the size of the atomic cloud, otherwise the use of def- 
inition (1351) becomes impractical due to roundoff errors. 
To summarize the discussion, the presented numerical re- 
sults only indicate the formation of the long-range corre- 
lations. The detailed description of the second-order cor- 
relation function in the limit of large separations can be 
obtained using a properly defined thermodynamic limit. 

V. SUMMARY 

The CM position measurement of trapped ultra-cold 
gases can become an important ingredient of various 
technologies based on ultra-cold atoms. The effect of 
such a measurement performed on quasi ID harmoni- 
cally trapped interacting Bose-gas has been studied in 
this article. The CM measurement back action disturbs 
the momentum, which due to the oscillations in the trap 
results in atomic delocalization. It has been shown that 



[1] J. Dunningham, K. Burnet, and W. D. Phillips, Phil. 
Trans. R. Soc. A 363, 2165 (2005). 

[2] T. Vaughan, P. Drummond, and G. Leuchs, Phys. Rev. 
A 75, 033617 (2007). 

[3] M. Antezza, L. P. Pitaevskii, and S. Stringari, Phys. Rev. 
A 70, 053619 (2004); D. M. Harber, et al, Phys. Rev. A 
72, 033610 (2005). 

[4] V. B. Braginsky and F. Ya. Khalili, Quantum measure- 
ment: Cambridge University Press, 1992. 

[5] F. Brennecke, et al, Nature 450, 268 (2007). 



the interaction-induced nonlinearity can to some extent 
resist this delocalization decreasing the rate of growth of 
the width of the average density profile. In other words, 
the atomic cloud size grows slowly for the stronger in- 
teracting gas. This result has first been obtained us- 
ing a semi-analytical mean-field approach with certain 
density-density type correlations ignored. Then the same 
conclusion has been obtained on a more rigorous basis 
performing numerical simulations based on positive P- 
representation of the many-atom density operator. In 
these simulations the atomic correlations have been taken 
into account, but the continuous density distribution has 
been approximated by a lattice model. 

Numerical simulations show that the averaged atom 
density profile or atomic cloud size remains almost un- 
changed for strongly interacting gas during the first quar- 
ter of the oscillation period. The numerical approach has 
allowed to calculate second-order correlation functions 
in the presence of the CM position measurement. The 
value of this correlation function for small atomic sepa- 
ration has been found to be greater than one for weakly 
and even strongly interacting gases, which is the atomic 
bunching. 

In addition, some preliminary conclusions can be made 
about the width of the instant atomic density profile. As 
discussed above, the average density profile of strongly 
interacting gas is almost unchanged in the beginning of 
the evolution. However, the CM position uncertainty 
being decoupled from the internal degrees of freedom 
constantly grows due to measurement back action noise. 
This is possible if the instant profile shrinks below its 
initial value. The latter conclusion is to be verified by 
direct calculations, which will be given elsewhere. 

VI. ACKNOWLEDGEMENTS 

The authors thank V.V. Kozlov for valuable discus- 
sions. One of the authors (DAI) is grateful to S. Wallen- 
towitz for the introduction to the theory of open many- 
atom systems. DAI also acknowledge Saint Petersburg 
government for the research grant. The authors acknowl- 
edge Saint Petersburg State University for a research 
grant. MSS acknowledges financial support from "Dy- 
nasty" foundation. 



[6] Y. Colombe, et al, Nature 450, 272 (2007). 
[7] S. Gupta, et al, Phys. Rev. Lett. 99, 213601 (2007). 
[8] W. W. Murch, et al, Nature Phys. 4, 561 (2008). 
[9] H. M. Wiseman and L. K. Thomsen, Phys. Rev. Lett. 
86, 12 (2001). 

[10] D. Ivanov and S. Wallentowitz, Phys. Rev. Lett. 93, 

260603 (2004). 
[11] F. Dalfovo, et al, Rev. Mod. Phys. 71, 463 (1999). 
[12] G. S. Ng, et al, New Journal of Physics 11, 073045 

(2009). 



9 



[13] H.-P. Breuer and F. Petruccione, The Theory of Open 

Quantum Sysetms: Oxford University Press, 2002. 
[14] P. D. Drummond and C. W. Gardiner, J. Phys. A: Math. 

Gen. 13, 2353 (1980). 
[15] A. Gilchrist, C. W. Gardiner, and P. D. Drummond, 

Phys. Rev. A 55, 3014 (1997). 
[16] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 
[17] I. B. Mekhov, C. Maschler, and H. Ritsch, Nature Phys. 

3, 319 (2007); I. B. Mekhov and H. Ritsch, Phys. Rev. 

Lett. 102, 020403 (2009). 
[18] Y. Colombe, et al, Nature 450, 272 (2007). 
[19] N. I. Kaliteevsky, Wave Optics: Vysshaya Shkola, 1978 

(in Russian). 



[20] C. M. Caves and G. J. Milburn, Phys. Rev. A 36, 5543 
(1987). 

[21] W. H. Press, et al., Numerical Recipes in C: The Art of 
Scientific Computing: Cambridge University Press, 1997. 

[22] P. Deuar and PD. Drummond, Phys. Rev. A 66, 033812 
(2002). 

[23] P.D. Drummond and I.K. Mortimer, J. Comp. Phys. 93, 
144 (1991). 

[24] M. Naraschewski, R.J. Glauber, Phys. Rev. A 59, 4595 
(1999). 

[25] Hereinafter we use h=l. 



